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The  Tsai  -  Wu  Constraint  is  Defined  in  ASTROS  as: 


Stiffness  Constraints 


37 


Sensitivity  Analysis 


c 

>  o 
■O  o 
c  c 


II  II 


D)  > 
ro  ro 


<D  ^ 
•  —  O) 

-Q - 

2  S 

0  Q 


u-  > 

ro 


CO  N. 
C  = 

o  05 
o  o 

(D  2a? 
-C  05 
■*"*  c 

o< 

W  o 
c  t; 

.2  Q. 

T3  c 

2  o 

O  O 


38 


Key  to  Performing  the  Approximate  Problem 
Computations  Can  Be  Intricate 


Architectural  Highlights 


39 


Provides  Programmer  with  Precise,  Explicit  Control 


Architectural  Highlights  -  Concluded 


2 

k_ 

o 

£ 

CL 

j 

TS 

'o 

* 

mmmm 

Cl) 

Q. 

C/D 

5 

1 

CL 

•  <0 
.9  E 

if)  k 
CO  tz 

0-0 

E_9> 

LU  < 


O 

<{>  * 

<D  c 

3  ;■§ 

■u  O 

o 


^  B 

o  w 

E  ~  .o 

I’cll 

^  LU 
—  « 

■—  Q)  .9 

pH 

2  jS  03 
5  d  w 

o  ro 

sla 

m'5g 

8  Is 

jcOC  < 


x  I — 

LU 

o  .52 


J2  .9 
o  c 
5  § 


^  E 

o  k 


—  to  to 
xi  c  cd  to 

CO  fl)  W  03 
to  0-5  2 
LU  O  O  LU 


to  T3 
CD  O 

SIN  ^ 


Relation  Of  Dynamic  Analysis  Sets 


co 

LL 

O 

Q 

*— » 

Q) 

to 

CO 

to 

>N 

(0 

co 

c 

< 

to 

4-* 

C 

o 

CD 

Q_ 

CO 

5 

CD 

X 

LU 

<0 

CD 

CO 

c 

■5 

o 

o 

Set 

O 

— 

15 

■o 

o 

2 

CD 

to 

c 

o 

Q. 

t  . 

to 

<D 

CD 

CO 

DC 

<*>* 

TJ 

o 

CD 

W_ 

b 

CD 

to 

C 

o 

Q. 

r  . 

to 

CD 

CD 

GC 

CO 

sz 

CO 

"D 

O 

2 
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Finite  Elements  -  The  Bar  Element 
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Non  -  Structural  Mass  is  Ignored 

Pin  Connection  and  Offset  Not  Supported 


Finite  Elements  -  The  Quadrilateral 
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Only  Isotropic  Materials  Supported 
Design  Variable  is  Element  Thickness 
Stress  Constraint  Based  on  the  Average  Shear  Stress 
at  the  Four  Nodes 
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Finite  Elements  -  The  Isoparametric 
Quadrilateral  Membrane  Element 
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Element  Warping  is  Allowed 
Anisotropic  Materials  Supported 
Design  Variable  is  Element  Thickness 

-  Comments  from  Triangular  Element  Apply 


Finite  Elements  -  The  Quadrilateral 
Shell  Element _ 
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Design  Variable  is  Element  Thickness 

-  Only  Membrane  Effects  Considered  in  Design 

-  Comments  from  Triangular  Element  Apply 


COMPARISON  OF  QUAD4  ElENENTS 
ASTROS,  NSC/NASTRAN1  AND  COSNIC  NASTRAN 
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UNIVERSAL  ANALYTICS,  INC 


Matrix  Assembly  -  Stage  One 
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Performed  as  a  Preface  Operation 


Matrix  Assembly  -  Stage  Two 
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Performed  as  a  Preface  Operation 

These  Matrices  are  Basic  to  ASTROS 
Sensitivity  Analysis 


Matrix  Assembly  -  Stage  Three 
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Static  Loads  -  Capability 
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{DPTH}.  =  X  Pjj  [Tee]j  {TGR|D  -  Tref }j 


Static  Loads  -  Assembly 
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Allows  for  Combining  of  Simple  Loads 
Aooommodafes  Design  Dependent  Loads 


Static  Analysis  -  Equations  Of  Motion 


6 


Where 


Static  Analysis  -  Solution  And 
Recovery 
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Recovery  to  g  -  set  is  Standard 


Static  Analysis  -  Strength  Constraint 
Evaluation 
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Stress /Strain  to  Global  Displacements 


Static  Analysis  -  Strength  Constraint 
Example 
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Static  Analysis  -  Strength  Constraint 
Sensitivity _ 
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Non  -  Zero  Only  f  Obtained  By 

For  Thickness  Constraints  )  \  Chain  Rule  Differentiation 


Gradient  Method  For  Sensitivity 
Analysis  -  Overview 
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Method  is  General 


Virtual  Load  Method  For  Sensitivity 
Analysis 
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Method  Not  Applicable  With  Inertia  on  Aerodynamic 
Terms 


Sensitivity  Analysis  -  Gradient  Method 


These  are  Pseudo  -  Load  Vectors  that  are  Designated  DP 


Sensitivity  Analysis  -  Gradient 
Method  (Concluded) _ 
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Modal  Analysis 
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Dynamic  Reduction 
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Dynamic  Reduction  -  Approximate 
Mode  Shapes 
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Starting  Vector 
Number  of  Iterates 
Shift  Point 

Rejection  of  Parallel  Vectors 
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Excellent,  But  Limited  Results 


Frequency  Constraint  Evaluation 
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Structural  Frequencies  Can  Be  Squeezed  Into 
a  Range  But  Not  Squeezed  Out 


Frequency  Constraint  Sensitivity 
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Steady  Aerodynamics 
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Aerodynamic  Influence  Coefficients 
Multiple  Mach  Numbers 


Steady  Aerodynamics 
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Rigid  Aerodynamic  Forces  are  Computed  For  a  Series  of 
Boundary  Conditions 

Aerodynamic  Influence  Coefficient  Matrix  Based  on  Linearized 
Pressure  Computation 


Unsteady  Aerodynamics 
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buDsomc  Aerodynamics 

No  Bodies  For  Supersonic  Aerodynamics 

Multiple  Lifting  Surfaces 


Unsteady  Aerodynamics 
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Design  Independent  Calculations  are  Performed  Only  Once 


Surface  Spline 
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Deformation  a  -  Aero  Model 

Slope  s  -  Structural  Model 
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Static  Aeroelastic  Analysis 
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A  New  Matrix  Combines  Structural  and  Aerodynamic  Stiffnesses 
[KAff]  =  [Kff  -  AlCSff] 


Static  Aeroelastic  Analysis 
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This  is  the  Basic  Equation  For  Static  Aeroelastic  Analysis 


Static  Aeroelastic  Analysis  - 


is  Determined  That  Provides  Required  Lift 
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Given  u2  and  5  ,  Recovery  of  Displacements  and 
is  Straightforward 
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Lift  Effectiveness  Constraint 
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upper  bouna,  Negative 
Can  Also  Be  Specified 
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Constraint  Form  and  Capabilities  Similar  to 
Lift  Effectiveness 


Static  Aeroelasticity  -  Sensitivity 
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Static  Aeroelasticity  -  Sensitivity 
Analysis 
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Configuration  Parameters  are  Invariant 
Acceleration  is  Invariant  and  Zero 
Changes  in  Displacement  are  Computed 


THE  EQUATION  REPRESENTS  A  SYNTHESIS  OF  NASTRAN  AND  FASTOP 
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GIVEN  p: 


Flutter  Analysis  -  Aerodynamic 
Interpolation 
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Quality  of  the  Interpolation  is  Assessed  By  the 
Procedure 


Flutter  Analysis  Algorithms 
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Flutter  Constraint  Properties 
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Gradients 
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Adjoint  Flutter  Vector  Utilized 

Similar  to  Frequency  Constraint  Sensitivity 

Procedure  is  Conceptually  Straightforward  But  Algebraically  Complex 


Automated  Design  -  Methods 
Of  Solution 
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computationally  ttticien 
Limited  in  Application 


Mathematical  Programming 
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Other  Algorithms  Could  Be  Readily  Substituted 


Reduction  In  The  Number  Of 
Desian  Variables 


200-300  Variables  Taxing  for  a  Micro  Computer 
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Inactive  Boundary  Conditions 
Inactive  Disciplines 
Inactive  Subcases 


The  Approximate  Design  Problem 
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Current  Value  of  the  Gradient  of  the  Active 
Constraints  with  Repect  to  the  Design 
Variables 


The  Approximate  Design  Problem 
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Redesign  Using  Inverse  Design 
Variables 
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Move  Limits  are  Imposed  on  the  Design  Variables  During  Redesign 


Redesign  Using  Direct  Design 
Variables 
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Termination  Criteria 
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Termination  Criteria  (Conducted) 
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Dynamic  Response  Analysis 
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Direct  and  Modal  Formulations  Available 
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Denotes  Direct  Input  Matrices 


Dynamic  Matrix  Assembly 
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Dynamic  Response  Solution 
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Transient  Response  with  a  Feedback  System 
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Damping  Matrix 


Response  as  a  Function  of  Gain 
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TIME  (SECS) 


Gust  Analysis 
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Existing  Code  was  a  Major  Resource  for  the  New 
Capability 
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The  Nuclear  Blast  Calculation 
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The  Aero  Preprocessor  For  Nuclear 
Blast  Calculation 
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Special  Treatment  is  Given  to  the  Early  Time  Loading  of  a  Box 
on  Itself 


The  Aero  Preprocessor  (Concluded) 
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Matrices  are  Converted  to  Generalized  Form  For  the  Response 
Calculation 


The  Blast  Response  Calculation 
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Vehicle  Deformations 
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User's  Interface  to  ASTROS 
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Executive  Control  Sequence  (MAPOL) 
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ASTROS  System  Generation 
Program,  SYSGEN 


SYSGEN  Output  As  User 
Documentation 
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Provides  Current  Listing  of  the  Standard  MAPOL 
Sequence 


Can  Be  Made  Available  "  On  -  Line 


ASTROS  User  Interface 


Similarities  Between  ASTROS  Input 
and  NASTRAN  Input 
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Selects  Output  Quantities 
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Transient 


Solution  Control  Contrasts  Between  NASTRAN  and  ASTROS 
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Solution  Control  Boundary  Condition  Type 
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Solution  Control  Boundary 
Condition  Definition 
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DAMPING  =  <  n  >  K2PP 

TFL  =  <  n  >  B2PP 

ESET  =  <n>  M2PP 


STATICS  Discipline  Options 


Each  Set  of  Options  Defines  a  Single  Load  Case  as  a 
Superposition  of  All  Load  Types 


STATICS  Discipline  Options 


co 

X S 

o  c 

CD  O 
CO  s 

He 
Z  o 

<° 

OC-D 

o-1 

o  © 

Q£ 


CL 

CO 

o 

z 

o 

o 

Q 


-L_  r- 

D)§ 

3  *■*= 
22 
SZ  O 
H  CO 


IdeT 

g-T3 

0)  a? 

©CO 

to  ^ 

II 
I  .a 
o| 

i5 

-ber 

COh- 

oz 

s° 

©D 

CO  ^ 
CO  C* 
0 


CO 


LU 


■e  fl 
jS  ©  — 

*=Q  O 
■£ 

O  D  O 
ZOO 


150 


No  STATICS  Options  are  Required  But  at  Least  One  of  MECH, 
GRAV  or  THERM  Must  Be  Present 


MODES  Discipline  Options 
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Flutter  Analysis  Requires  That  the  Eigenvalue  Extraction  Method 
Be  Specified  in  the  Boundary  Definition 


TRANSIENT  Discipline  Options 
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FREQUENCY  Discipline  Options 
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MODAL  Analyses  Require  BOUNDARY  METHOD  =  <  n  >  to 
Perform  Real  Eigenanalysis 


BLAST  Discipline  Options 
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MODAL  Analyses  Require  BOUNDARY  METHOD  =  <  n  >  for 
Real  Eigenanalysis 


Combining  Disciplines  In  A  Single 
Boundary  Condition 
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Precluded  R  -  Allowed  With  Restrictions 


Multidisciplinary  Restrictions 
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Only  One  Modal  Analysis  Allowed  in  a  Boundary  Condition 

TRANSIENT,  FREQUENCY  and  BLAST  are  Limited  to  One 
Analysis  Each  Per  Boundary  Condition 
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Once  Selected,  an  Output  Request  Remains  in  Force  at  or 
Below  that  Level  in  the  Hierarchy  Until  Overridden 


Solution  Control  Output  Requests 
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POLAR  Selects  Magnitude  /  Phase  of  Complex  Quantities 
<  form  >  is  Ignored  for  Real  Quantities 


Solution  Control  Output  Requests: 
Response  Quantity  Options 
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Solution  Control  Output  Requests: 
Response  Quality  Options  (Concl'd) 
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Discipline  Independent 


Solution  Control  Output  Requests: 
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That  the  Print  Request  Applies  to  NO 
Subcases 


Solution  Control  Output  Requests: 
Common  Pitfall 
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This  Example  Results  in  DISP  =  5  for  Both 
STATICS  and  MODES 
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of  the  1st  Discipline 


Solution  Control  Example 
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TRANSIENT  MODAL  (DLOAD  =  500,  TSTEP  =  501) 
PRINT  DISPLACEMENT  =  ALL,  TIME  10,  STRAIN  =  12 


Bulk  Data  Packet 


170 


ENDDATA 


Bulk  Data  Entry  Formats 
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Bulk  Data  Entries: 
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MKAER01,  ,  0.3,  0.5 . ABC  ,  +BC,  0.01,  0.05,  0.1,  0.2 

MKAEROI,  ,  0.3,  0.5 . 0.01,  0.05,  0.1,  0.2 


Bulk  Data  Entries:  Data  Fields 
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Blank  Fields  are  Automatically  Replaced  with  0,  0.0  or 
Depending  on  the  Reid  Type 


ASTROS/NASTRAN  Bulk  Data 
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Fields 


Bulk  Data  Entry  Differences 
Between  ASTROS  And  NASTRAN 
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-  Design  Constraints 

-  Steady  Aerodynamics  Model  Geometry 

-  Discipline  Data  for  New  Disciplines 

Remainder  are  Changed  to  Some  Degree 


ASTROS  Modifications  To  Existing 
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Multidisciplinary  Analysis  Changes 


Shape  Function  Desian  Variable 
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The  Design  Variable  in  ASTROS 


Design  Variables 


Unique  Physical  Design  Variable 
Linkina 


83 


ONE  ELEMENT/LAYER  FOR  EACH  GLOBAL  VARIABLE 


Shape  Function  Design  Variable 
Linking 


o 


o 


w 


05 


00 


CO 


in 


co 


cvi 


85 


MULTIPLE  ELEMENTS  FOR  EACH  GLOBAL  VARIABLE  AND/OR 
MULTIPLE  VARIABLES  FOR  EACH  ELEMENT/LAYER 
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The  DCONTHK  Bulk  Data  Entry  Has  Been  Defined  for  this 
Purpose 
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Shape  Function  Linking  Enforces  an  Initially  Uniform  Local 
Property  Distribution 


Limitations  In  Design  Variable 
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Two  Layers  of  a  Composite  Element  Cannot  Be  Linked  to 
the  Same  Design  Variable 


Design  Constraints  In  ASTROS 
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Symmetric  Analyses  Only 
Antisymmetric  Analyses  Only 


Stress/Strain  Constraint  Specification 
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VMISES  for  von  Mises  Stress  Criterion 
TSAIWU  for  Tsai  -  Wu  Strength  Ratio 
STRAIN  for  Principal  Strain  Constraint 


DCONSTR  And  The  MAT1 
Material  Property 
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DCONSTR  And  The  MAT2 
Material  Property 
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DCONSTR  And  The  MAT8 
Material  Property 


I- 


o 

05 

00 

1- 

Z 

o 

o 

O 

CD 

< 

+ 

K 

Z 

o 

o 

Li. 

HI 

D 

+ 

o 

X 

X 

CO 

in 

o 

o 

w 

c/5 

CO 

+ 

N 

cvi 

O 

CO 

+ 

in 

o 

>- 

CM 

+ 

00 

N 

5 

CO 

+ 

CO 

£ 

CM 

+ 

cvi 

co 

CM 

5 

CO 

+ 

cvi 

£ 

M- 

+ 

in 

CM 

T— 

D 

Z 

CO 

o 

X 

TT 

+ 

1 

■m- 

CO 

CM 

CM 

LU 

CO 

+ 

LL 

UJ 

X 

H 

o 

in 

m 

i — 

LLJ 

CO 

+ 

o 

CO 

a 

CO 

I 

in 

T“ 

CM 

T— 

U. 

Q 

2 

g 

< 

CO 

1 

00 

CM 

0 

E 

00 

s 

2 

00 

< 

2 

H 

Z 

o 

o 

O 

CO 

+ 

H- 

Z 

o 

o 

Li. 

LU 

O 

+ 

1 

i 


£ 


I 

1 

I 


S 


O 

1 

O  O 

fl  < 

I  8 

I  g 

w 

V) 

W  fl) 
2?  £ 
tn 

07 

.Q 

iZ 


c 

o 


a> 

> 

CO 

c 

CO 


c  c 
o  o 

'8  ’8  5 

£  £  -8 

a  a  | 

E  E  o 
o  o  = 
O  O  < 

X3  O  S 

S  S  S 


c  CO 
O 

M  tt  W 

C  C  <D 

0  CD  JC 

h  h  W 


o  o 
X  > 

-  —  co  £! 
X  >■  co  E 


194 


For  Principal  Strain 

Xt ,  Xc  -  Tension  and  Compression  Strain  Allowables  in  Microunits/Unit 

Y , ,  Y„  .SS.F12  -  Not  Used 


Additional  Stress/Strain 
Constraint  Information 
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An  Element  Need  Not  Be  Designed  to  Be  Constrained 


Displacement  Constraint  Specification 


LOWER  'A  y/s  UPPER 


Frequency  Constraint  Specification 
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Cannot  Use  Multiple  Constraints  to  Exclude  a 
Frequency  From  a  Region 


Flutter  Constraint  Form 
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Aileron  Effectiveness  Constraint 
Specification 


zoo 


Required  Effectiveness  May  Be  Positive 
or  Negative  so  That  a  Reversed  Aileron 
Condition  May  Be  Imposed 


Lift  Effectiveness  Constraint 
Specification _ 
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Positive  or  Negative 


Unsteady  Aerodynamic  Parameters 
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Velocity  Field  is  Redundant 
Symmetry  Fields  Limit  Multidisciplinary 
Capability 


Unsteady  Aerodynamic  Lifting  Surfaces 
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Panel  Boundaries  -  Edges  are  Parallel  to  the  Flow 


Unsteady  Aerodynamic  Body  Connection 


Leading  Point  of  Body 


Interference 
Elements  Us 


Flutter  Analysis  Conditions 


Multidisciplinary  Analysis  Requires  Symmetry  Condition 
as  Part  of  the  Specification 


Bulk  Data  Entries  For  Steady 
Aerodynamics 


<D 

■o 

> 

O 

CL 

03 

To 

Q 


°  £ 
O  O 


207 


Paneling  Data  Define  the  Mathematical  Representation  in 
USSAERO 


Steady  Aerodynamic  Parameters 


*—  Reference  Coordinate  System 

NOT  OPERATIONAL 

Aerodynamic  Coordinate  System  IN  ASTROS 


ID  For  Chordwise  Paneling  Cuts 
Chord  Length 

Edge  Location 


Airfoil  Properties  (Concluded) 
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Lower  =  -100  z.  /c  Camber 


Steady  Aerodynamic  Lifting  Surface 
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Chordwise  Cuts  are  in  Percent  Chord 
Spanwise  Cuts  are  in  Physical  Coordinates 


Body  Properties 


Station  Location 


Steady  Aerodynamic  Body  Surface 
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Limits  on  Paneling  Data  in 
USSAERO 
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Control  Surface  Definition 
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Not  Used 


Flight  Condition  Parameters 


TRIM  ID  From  Solution  Control 
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The  ASTROS  Input  File  Processor 
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Must  Be  Called  By  the  Executive  for  Every 
Execution 


The  ASTROS  Bulk  Data  Template 
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User  Output  From  ASTROS 
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System  Controlled  Output  From 
ASTROS 
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ystem  and  user  Error  Messages 
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Termination  States  and  Timing  Summary 


Boundary  Condition  Default 
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Active  Constraint  Summa 
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— 1.45733E- 
6.48499E- 


Approximate  Optimization  Summary 
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Design  Information:  Iteration 
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THE  FINAL  OBJECTIVE  FUNCTION  VALUE  IS: 

FIXED  =  2 . 90300E+01 

+  DESIGNED  =  2.64772E+01 


Final  Design  Information:  Design 
Variables  _ 


9  UNIQUE  PHYSICAL  1 . 80388851E+01  1.000E-01  1.000E+04 

10  UNIQUE  PHYSICAL  6 . 99411273E-01  1.000E-01  1.000E+04 


Order  Of  Output  For  Selected  Quantities 
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Within  Each  Quantity,  The  Disciplines  are  Treated: 

(1)  STATICS  (5)  TRANSIENT 

(2)  MODES  (6)  FREQUENCY 

(3)  SAERO  (7)  BLAST 

(4)  FLUTTER 


OFP  Example  -  Stability  Derivatives 


OFP  Example  -  Stress  Output 
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em  And  User  Error  Messages 
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(2)  User  Information 

(3)  User  Warning 

(4)  User  Fatal 


Determining  The  Source  Of  An  Error 
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0  ***FREE****  -1  769850  769844 

769850  RELINDEX  0  770368  512  DBB1 

770368  PGMTST  769850  770886  512  DBB] 

770886  RELSCHEM  770368  772940  2048  DBB1 


Executive  Controlled  Output 
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General  Relational  Print  Utility 
General  Unstructured  Entity  Print  Utility 


Optional  Print  Arguments  For  IFP 
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Echo  Only  to  Punch  File 

Echo  to  Both  Output  and  Punch  Files 


Optional  Print  Argument  For  DESIGN 


•rder  of  Constraints  May  Not  Match  Other  ASTROS  Output 
iternal  Scaling  Further  Modifies  the  Output 


ASTROS  Executive  System 
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10-15-87-05-EJ 


TheMAPOL  Packet 


MAPOL  Edit  Commands 


Structure  Of  The  Standard  MAPOL  Sequence 
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Features  Of  The  MAPOL  Language 
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Features  Of  The  MAPOL  Language 
(Concluded)  _ _ 
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TRNSPOSE 


User  Supplied  MAPOL  Program  Example 
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Modifying  The  Standard  MAPOL 
Sequence 


=  •§ 

8 1 

O  CO 


-Q  A 


£  o> 
a>  c 

cc  ~ 
(/> 


CO  (0  CO 
LL) 

(_  LU  O 

f  h  < 

5  y  i 

^2  Dj  lu 
_  O  CC 


C  3 
(0  Q. 


a>  z 
£1  111 


E  O  8 
3«C 


a>  §" 
£  co 


o  W 

<  g 

O  <0 


245 


Editing  Must  Be  Done  in  Increasing  Line  Number 
Order 


ASTROS  Executive  Sequence  - 
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Typical  Replacements  to  Standard  Sequence  Involve  : 

-  Restart  to  Compute  and  Print  Additional  Data 

-  Special  Purpose  Analyses 
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UNIVERSAL  ANALYTICS,  INC. 
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UNIVERSAL  ANALYTICS,  INC 


ASTROS  SOFTWARE  DESIGN 
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UNIVERSAL  ANALYTICS,  INC 


UNIVERSAL.  ANALYTICS,  INC 


ASTROS  SOFTUARE  DESIGN 


MMVTICt,  INC 


ASTROS  SOFTWARE  DESIGN 
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Using  Arithmetic  Expressions 
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The  JOIN  Operation 
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The  SELECT  COLUMNS  Command 


cq 

I 

o 

o 

in 

•*  cq 

0  2 

o  o 

2  2 

2  in 

H 

H 

o 

2 

2 

2 

2  2 

o 

O 

O 

01 

05 

CN  H 

*  CQ 

Eh 

EH 

H  H 

< 

< 

0 

—  w 

CQ 

CQ  O 

0  2 

2 

2  O 

2  04 

H 

H  O 

H 

0 

0  O 

2  w 

W 

w  o 

Eh  3 

CQ 

CQ  • 

(0  CQ 

D 

1 

1  ** 

CQ  O 

2  Q 

H 

H 

£ 

O 

O 

5  q 

0  O 

0  O 

2  O 

2  O 

o  a 

H  O 

H  O 

u  2 

0$  O 

2  o 

Eh  • 

Eh  • 

fcj  u 

CQ  CN 

CQ  co 

W  0 

3  2 

H 

cn 

w 

o 

o 

CQ  X 

2  o 

2  O 

H 

2  o 

£  o 

A  2 

5  o 

R  o 

w  Eh 

3  o 

3  o 

0 

o  • 

o  • 

H  X 

0  H 

0  CM 

CQ 

l 

i 

l 

i 

l 

l 

i 

l 

g 

i 

l 

_ . 

i 

IN 

o 

i 

o 

’ — ' 

0 

i 

o 

l 

o 

2 

in 

i 

o 

H 

i 

I 

o 

CQ 

cq" 

1 

i 

o 

2 

2 

i 

O 

O 

l 

2 

2 

l 

1 

o 

w 

in 

l 

o 

2 

i 

o 

w 

i 

o 

2 

i 

o 

O 

i 

• 

H 

1 

o 

0 

CQ 

1 

0 

H 

i 

2 

0 

i 

l 

o 

2 

2 

l 

o 

O 

2 

l 

o 

2 

i 

o 

2 

3 

i 

1 

o 

* 

cq 

1 

l 

O 

i 

CQ 

2 

i 

l 

l 

o 

5 

i 

o 

rij 

l 

rH 

o 

o 

(jj 

l 

o 

o 

2 

l 

g 

o 

l 

p»j 

• 

EH 

*. 

l 

5 

CO 

0 

0 

l 

0 

0 

l 

o 

3 

2 

l 

0 

w 

l 

o 

CQ 

X 

i 

o 

H 

l 

CN 

o 

A 

2 

l 

o 

W 

» 

2 

o 

0 

< 

1 

o 

• 

H 

2 

1 

2 

CN 

335 


UNIVERSAL  ANALYTICS,  INC. 


336 


UNIVERSAL  ANALYTICS,  INC. 


The  CREATE 


337 


UNIVERSAL  ANALYTICS,  INC. 


338 


UNIVERSAL  ANALYTICS,  INC. 


The  INSERT  Commands 
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Removinq  Data  from  CADDB 
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COLUMN  Command 
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Rectangular  Wing  Box  Model 
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Design  Requirements  For  Intermediate  Complexity  Wing 
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Design  Cases  For  Intermediate  Complexity  Wing 
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utter  Contraint  is  Impose 
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ICW  Strength  &  Flutter  Design 

Ply  Counts  for  the  -45°  Laminate 
Three  Numbers  are : 
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Structural  Model  of  the  N372-4  Fighter  Wing 
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1,494  Constrained 
1,506  Unconstrained 

3,000  Total 
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N 372-4  Example  —  Iteration  History 


Iteration  Number 


ASTROS  User  Training  Workshop 

Problem  Set  Definitions 

Problem  Set  #1:  Space  Truss 

1-1 :  Modal  Analysis  using  Guyan  and  Generalized  Dynamic  Reduction 

1- 2:  Optimization  for  first  two  Natural  Frequencies 

Problem  Set  #2:  Rectangular  Wing 

2- 1 :  Static  Analysis  for  Tip  Load 

2-2:  Optimization  for  Stress  constraints 

2-2:  Optimization  for  Stress  constraints  with  Shape  Functions 

2-4:  Static  Analysis  with  Inertia  Relief 

2-5:  Aeroelastic  Trim  for  Wing  in  Straight,  Level  Flight 

2-6:  Aeroelastic  Trim  for  Wing-Tail  combination  in  Pull-up  maneuver 

2-7 :  Aeroelastic  Analysis  for  Roll  maneuver 

2-8:  Optimization  for  Stress  and  Tip  Twist 

2-9:  Optimization  for  Stress,  Tip  Twist,  and  Lift  Effectiveness 

2- 10:  Optimization  for  Stress,  Tip  Twist,  Lift  and  Aileron  Effectiveness 

Problem  Set  #3:  Cantilvered  Plate 

3- 1 :  Static  and  Modal  Analyses 

3-2:  Transient  Analysis 

3- 3:  Subsonic  and  Supersonic  Flutter  Analysis 

Problem  Set  #4:  Swept  Wing 

4- 1 :  Static  Analysis  for  Gravity  Load 

4-2:  Modal  Analysis 

4-3:  Optimization  for  Stress  and  Frequency  constraints 
4-4:  Subsonic  Flutter  Analysis 

4- 5:  Supersonic  Air  Loads 

Problem  Set  #5:  Plane  Frame 

5- 1 :  Optimization  of  40  member  Plane  Frame  for  Stress  and  Displacement  Constraints 


Workshop  Requirement  is  to  complete  10  of  the  above  21  problems, 
including  at  least  one  from  each  problem  set. 
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ASTROS  User  Training  Workshop 
Problem  Set  #1 :  AC  OSS  Space  Truss 


The  Active  Control  Of  Space  Structures  (ACOSS)  model  II  was  developed  by  the  Charles  Stark 
Draper  Laboratory.  The  structure  consists  of  two  subsystems:  (1)  the  optical  support  structure 
and  (2)  the  equipment  section.  The  two  are  connected  by  springs  at  three  points  to  allow 
vibration  isolation  (Figure  1-1).  For  this  problem  set  disregard  the  equipment  section  at  the  base 
and  consider  the  optical  support  structure  fixed  at  the  three  connection  points.  The  finite  element 
model  for  this  modified  ACOSS  II  (Figure  1-2)  has  33  nodes  (90  degrees  of  freedom),  18 
concentrated  masses,  and  113  rod  elements  made  of  graphite  epoxy  given  (Table  1-1)  with  initial 
areas  of  10.0  in2  for  the  truss  members.  The  grid  points,  masses,  and  element  connectivities  are 
given  in  attachment  1. 

For  this  initial  design  the  first  three  frequencies  of  the  modified  ACOSS  II  truss  are:  1.21,  2.71, 
and  4.21  hz. 

1-1)  Verify  the  first  three  frequencies.  Compare  Guyan  Reduction  (omit  degrees  of  freedom 
for  nodes  without  concentrated  masses)  to  Generalized  Dynamic  Reduction  for 
reducing  the  size  of  the  problem  before  applying  Givens  method. 

1-2)  Design  the  truss  for  minimum  weight  while  raising  the  fundamental  frequency  to  2  hz 
and  maintaining  at  least  a  1  hz  separation  of  the  fundamental  mode  from  the  remaining 
modes.  Use  a  minimum  gage  size  of  0.01  in2  for  the  truss  elements.  What  are  the  first 
three  frequencies  and  weight  for  the  final  design? 


Young’s  Modulus 

18.5  x  106  psi 

Weight  Density 

0.055  lb/in3 

Table  1-1:  Material  Properties  for  Epoxy 


ACOSS  Space  Truss 
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ACOSS  Space  Truss 
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Problem  Set  #1 


•  Lumped  mass  location 
- Added  support  rods 


Figure  1-2:  Finite  Element  Model  for  Modified  ACOSS  II 

ACOSS  Space  Truss  Problem  Set  #1 
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ASTROS  User  Training  Workshop 
Problem  Set  #2:  Rectangular  Wing 


Structural  Wing  Box  Model 

A  simple  three-spar  rectangular  wing  box  is  shown  in  Figure  2-1.  The  semi-span  is  60”  and  the 
chord  of  the  structural  box  (distance  from  front  to  rear  spar)  is  20”.  A  200  lb  concentrated  mass 
with  a  moment  of  inertia  about  the  span  axis  of  22,500  lb-in 2  at  the  root  of  the  front  spar 
represents  the  fuselage  mass.  The  wing  is  made  of  aluminum  (material  properties  given  in  Table 
2-1).  The  structural  box  is  modeled  using  quadrilateral  membrane  elements  for  the  cover  skins, 
shear  panels  for  the  spars  and  ribs,  and  rod  elements  for  the  vertical  posts  (cross-sectional 
properties  given  in  Table  2-2). 


K. 

Figure  2-1:  Rectangular  Wing  Structural  Box 


Young’s  Modulus 

10.0  x  106  psi 

Poisson’s  Ratio 

0.3 

Weight  Density 

0.1  Ib/in3 

Tensile  Strength 

20.0  ksi 

Compressive  Strength 

15.0  ksi 

Shear  Strength 

12.0  ksi 

Table  2-1:  Material  Properties  for  Aluminum 


Rectangular  Wing 
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Membrane  Thicknesses 

0.20  in 

Shear  Thicknesses 

0.05  in 

Rod  Areas 

0.01  in2 

Table  2-2:  Cross-Sectional  Properties 

2-1)  Perform  a  static  analysis  of  the  wing  when  it  is  subjected  to  100  lb  load  applied 

vertically  at  each  of  the  six  grid  points  at  the  wing  tip.  Consider  the  wing  cantilevered 
at  the  root.  Find  the  displacements  of  each  grid  point  and  the  stresses  in  the  cover  skins. 

Design  Model 

Use  design  variable  linking  to  define  four  design  variables  that  control  the  thicknesses  of  the 
eight  cover  skin  membrane  elements.  Each  design  variable  controls  a  group  of  two  membrane 
elements,  one  fore  and  one  aft,  so  the  design  variables  control  the  (1)  outboard  upper  skin,  (2) 
outboard  lower  skin,  (3)  inboard  upper  skin,  and  (4)  inboard  lower  skin.  The  spars,  ribs,  and 
posts  remain  Fixed  (are  not  designed). 

2-2)  Optimize  the  structural  weight  of  the  cover  skins  subject  to  stress  constraints  (Table  2- 
1)  on  the  cover  skins  only  .  Find  the  optimum  weight  and  design  variable  values  for  the 
boundary  conidtion  and  static  mechanical  load  given  in  problem  2-1. 

2-3)  Repeat  problem  2-2  using  shape  function  design  variable  linking.  Use  a  constant 
thickness  (initially  0.1”)  and  a  spanwise  linear  shape  (initially  0.075”  inboard  and 
0.025”  outboard)  for  each  of  the  upper  and  lower  skins  (4  design  variables). 

2-4)  Perform  a  static  analysis  of  the  wing  with  inertia  relief  for  the  static  load  given  in 
problem  2-1.  Use  multipoint  constraints  to  rigidly  connect  the  six  grid  points  at  the 
wing  root  to  the  root  of  the  center  spar  midway  between  the  top  and  bottom  skins.  The 
fuselage  mass  is  associated  with  this  “aerodynamic  reference  point”  with  an  offset  to 
locate  it  at  the  mid-surface  of  the  front  spar.  Support  the  aerodynamic  reference  point 
in  vertical  translation  (plunge)  and  find  the  displacements  and  accelerations. 

Steady  Aerodynamic  Panel  Model — Wing 

The  aerodynamic  planform  (Figure  2-2)  for  the  wing  has  a  30”  chord  and  60”  semi-span.  The 
structural  box’s  front  and  rear  spar  are  located  at  the  13.33%  and  80%  chord  locations, 
respectively.  The  aerodynamic  box  pattern  shown  in  Figure  2-2  has  four  chordwise  cuts  at  0%, 
20%,  80%,  and  100%  of  the  chord  and  five  equal  spanwise  cuts.  The  airfoil  shape  given  in  Table 
2-3  is  for  a  symmetric  airfloil  (no  camber)  with  a  leading  edge  radius  of  1.667%g.  An  aileron  is 
defined  by  the  two  outboard  trailing  edge  boxes  of  the  wing. 

2-5)  Find  the  trimmed  angle  of  attack,  displacements,  and  accelerations  for  symmetric  level 
flight  (lg  load  factor)  at  Mach  0.8  and  a  dynamic  pressure  of  6.5  psi.  Use  the  wing  only 
and  spline  the  aerodynamic  boxes  to  the  upper  sufurface  structural  grid  points. 

Compare  the  lift  coefficient  to  the  theoretical  value  for  a  thin  airfoil  wing. 


Rectangular  Wing 
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Figure  2-2:  Aerodynamic  Planform  for  Rectangular  Wing 


Rectangular  Wing 
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Chord  station 

Half  thickness 

0. 

0. 

10. 

1.667 

25. 

1.667 

50. 

1.667 

80. 

1.667 

100. 

0. 

Table  2-3:  Airfoil  Shape  (all  units  in  %  chord) 

Steady  Aerodynamic  Panel  Model — Tail 

The  trailing  edge  of  the  elevator  in  figure  1  is  80”  from  the  aerodynamic  reference  point  (center 
spar).  It’s  root  chord  is  15”,  its  tip  chord  is  10”,  and  its  semi-span  is  20”.  The  aerodynamic  box 
pattern  for  the  horizontal  tail  shown  in  Figure  2-2  (3  equal  spanwise  cuts,  and  4  chordwise  cuts 
at  0%,  20%,  80%,  and  100%  of  the  chord).  The  thickness  distribution  is  the  same  as  for  the  wing 
as  given  in  Table  2  with  a  leading  edge  radius  of  1%  An  elevator  is  defined  by  the  two  trailing 
edge  boxes  on  the  tail. 

2-6)  For  the  wing-tail  combination  find  the  trimmed  angle  of  attack,  elevator  deflection,  and 
tip  displacement  for  a  symmetric  8g  pull-up  maneuver  defined  by  the  flight  condition  in 
Table  2-4.  Trim  the  aircraft  for  lift  and  pitching  moment.  Support  the  structural  model 
for  pitch  and  plunge  rigid  body  modes. 


Mach  number 

0.8 

Dynamic  Pressure 

6.5  psi 

Pitch  Rate 

15.7  deg/sec 

Velocity 

487.4  knots 

Table  2-4:  Flight  Condition 

2-7)  Find  the  rigid  and  flexible  stability  derivatives  for  an  anti-symmetric  roll  maneuver 
using  the  flight  condition  in  Table  2-4. 

For  the  following  problems  optimize  the  structure  using  the  design  variables  for  problem  2-2. 
For  the  final  design  in  problems  2-8  through  2-10  find  the  values  of  constraints  not  imposed 
during  that  optimization.  Which  constraint(s)  are  critical  in  driving  the  design? 

2-8)  Optimize  the  structure  for  the  symmetric  8g  pull-up  of  problem  2-6.  Impose  stress 
constraints  (Table  2-1)  on  the  skins  and  a  maximum  tip  rotation  of  1  degree.  Subtract 
the  rotation  of  the  support  point  from  the  relative  rotation  of  the  tip  to  calculate  the  pure 
elastic  twist  of  the  tip. 

2-9)  Repeat  problem  2-8  with  the  addition  of  maximum  lift  effectiveness  of  1 .60. 

2-10)  Repeat  problem  2-9  with  the  addition  of  a  minimum  aileron  effectiveness  of  0.30  for 
the  anti-symmetric  roll  maneuver  of  problem  2-7. 
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Problem  Set  #3:  Cantilevered  Plate 


Structural  Plate  Mode! 

The  cantilevered  aluminum  plate  in  Figure  3-1  is  a  parallelogram  (constant  chord,  no  taper) 
swept  15°  with  a  uniform  thickness  of  0.041”.  The  tip  is  5.52”  from  the  cantilevered  root  and  the 
unswept  width  (chord)  is  2”.  The  material  properties  are  given  in  Table  3-1.  The  finite  element 
model  consists  of  a  course  3x5  mesh. 


Figure  3-1: 

Cantilivered  Plate 

Young’s  Modulus 

10.0  x  106  psi 

Poisson’s  Ratio 

0.33 

Mass  Density 

2.59  x  104  lb-sec2/in4 

Table  1 — Material  Properties  for  Magnesium 

3-1)  Perform  a  static  analysis  to  determine  the  displacements  for  a  unit  moment  applied  at 
the  one-third  chord  location  of  the  tip.  Next  find  the  first  three  natural  modes. 

3-2)  Determine  the  transient  response  in  the  time  domain  for  the  time-dependent  load  given 
in  Figure  3-2  applied  at  the  two  free  comers  of  the  plate  in  the  transverse  direction. 
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Time  Dependent  Load  Applied  to  Plate 


Time  (sec) 

Figure  3-2:  Time  Dependent  Load  for  Cantilvered  Plate 
Unsteady  Aerodynamic  Model 

An  unsteady  aerodyanmic  model  consists  of  5  span  wise  and  10  chordwise  aerodynamic  boxes, 
equally  distributed  over  the  planform  defined  by  Figure  3-1. 

3-3)  Determine  the  flutter  speed  for  the  subsonic  and  supersonic  flight  conditions  defined  in 
Table  3-2. 


Mach  number 

0.45 

3.0 

Air  Density  Ratio 

0.9676 

0.3913 

Reference  Density 

11.46  x  10‘6  lb-sec2/in4 

11.46x  10'6  lb-sec2/in4 

Table  3-2 — Flight  Condition 
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Problem  Set  #4:  Swept  Wing 


Structural  Wing  Box  Model 

The  planform  for  a  swept  wing  in  Figure  4-1  shows  the  top  skin  of  a  structural  model.  The 
structural  model  divides  the  structural  box  into  six  equally  spaced  spanwise  bays  and  two  equal 
chordwise  segments.  The  skins  on  both  the  upper  and  lower  surface  are  modeled  as 
isoparametric  quadrilateral  membrane  elements.  The  ribs  and  spars  are  modeled  using  shear 
panels  with  rod  elements  for  the  spar  caps.  Rod  elements  are  also  used  as  posts  connecting  all 
upper  and  lower  surface  nodes.  This  results  in  57  rod  elements,  24  quadrilateral  membrane 
elements,  and  32  shear  panels.  The  material  properties  are  given  in  Table  4-1,  cross-sectional 
properties  in  Table  4-2.  The  six  nodes  at  the  wing  root  are  fixed  (cantilvered). 

4-1)  Perform  a  statics  analysis  for  a  4g  vertical  gravity  load  and  find  the  displacements  and 
stresses. 

4-2)  Perform  a  modal  analysis  to  determine  the  first  five  normal  modes  of  the  structure.. 

Design  Model 

The  design  model  consists  the  sizes  of  the  skins,  spar  webs,  spar  caps,  and  wing  ribs.  Use  design 
variable  linking  to  couple  elements  in  each  of  three  spanwise  segments  of  the  structure  (12 
design  variables). 

4-3)  Optimize  the  structural  box  (excluding  posts)  subject  to  stress  constraints  on  the  wing 
skins  (24  constraints)  for  the  4g  gravity  load  and  a  1 .5  hz  lower  bound  frequency 
constraint  on  the  first  bending  mode. 
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Young’s  Modulus 

10.0  x  106  psi 

Poisson’s  Ratio 

0.3 

Weight  Density 

0.1  lb/in3 

Tensile  Strength 

60.0  ksi 

Compressive  Strength 

50.0  ksi 

Shear  Strength 

30.0  ksi 

Skin  Thicknesses 

0.16  in 

Rib  Shear  Thicknesses 

0.16  in 

Spar  Shear  Thicknesses 

0.32  in 

Post  Rod  Areas 

0.3  in^ 

Spar  Cap  Rod  Areas 

2.0  in^ 

Table  4-2:  Cross-Sectional  Properties 


Aerodynamic  Models 

The  planform  for  aerodynamic  and  structural  models  are  shown  together  in  Figure  4-2.  Both  the 
steady  and  unsteady  aerodynamic  models  represent  the  wing  as  a  flat  plate  with  50  boxes  per 
surface.  The  unsteady  model  has  ten  equally  spaced  spanwise  boxes  and  five  chordwise  boxes, 
while  the  steady  model  has  its  chordwise  boxes  spaced  in  a  cosine  distribution  (xj  =  C[  1  - 
cos(ip/5)]/2).  The  steady  model  has  a  horizontal  stabilizer  to  enable  trim  for  both  lift  and 
pitching  moment.  Like  the  wing,  the  tail  is  represented  as  a  flat  plate  with  ten  equally  spaced 
spanwise  boxes  and  five  chordwise  boxes  distributed  using  a  cosine  distribution.  The  last  two 
boxes  in  each  chordwise  strip  are  used  to  represent  an  elevator.  No  structure  is  associated  with 
this  tail  panel.  Both  aerodynamic  wing  models  transfer  the  forces  to  the  structural  nodes  on  the 
upper  surface  of  the  structural  box  with  a  linear  surface  spline.  The  tail  forces  for  the  steady 
aerodynamic  model  are  rigidly  transferred  to  the  center  root  of  the  structural  box. 

4-4)  For  flight  condition  1  in  Table  4—3  (Mach  0.8  at  sea  level)  and  the  strucutral  design 
point  given  in  Table  4-2  determine  whether  the  wing  flutters. 

4-5)  For  flight  condition  2  in  Table  4-3  (Mach  1.25  at  25,000  feet),  find  the  trimmed  angle 
of  attack  and  elevator  deflection  for  a  symmetric  4g  pull-up  maneuver. 


Flight  Condition 

1  (unsteady) 

2  (steady) 

Mach  number 

0.8 

1.25 

Load  Factor 

l.Og 

4.0g 

Elevation 

0.  ft 

25.0  x  103  ft 

Air  Density  Ratio 

1.0 

0.4486 

Dynamic  Pressure 

5.959  psi 

Velocity 

530.0  knots 

752.6  knots 

Pitch  Rate 

4.354°/sec 

Table  4-3:  Flight  Conditions 
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Problem  Set  #5:  Forty  Member  Plane  Frame  Optimization 

The  forty  member  plane  frame  shown  in  Figure  5-1  is  subjected  to  the  three  independent  loading 
conditions  described  below: 

Loading  condition  1: 

Distributed  vertical  load  O'  direction)  of  2.4  kips/ft  on  the  intermediate  levels  and  1.4  kips/ft  at 
the  top  level.  These  distributed  loads  are  approximated  by  nodal  forces  (25%  of  the  distributed 
loads  at  the  end  points  and  50%  at  the  mid  span). 

Loading  condition  2: 

Horizontal  forces  from  the  left  as  shown  in  figure  5-1  (solid  arrows)  plus  75%  of  loading 
condition  1. 

Loading  condition  3: 

Horizontal  forces  from  the  left  as  shown  in  figure  5-1  (dashed  arrows)  plus  75%  of  loading 
condition  1. 

Using  design  variable  linking  define  20  design  variables  in  which  the  two  horizontal  beams  at 
each  level  are  grouped  into  one  variable  and  similarly  the  the  two  vertical  members  at  each  level 
are  grouped  into  one  variable.  The  material  properties  and  the  intial  design  variables  are  given  in 
Table  5-1. 

Design  Problem:  Optimize  the  structure  for  minimum  weight  subjected  to  horizontal 
displacement  constraints  of  +2”  and  -2  inches  at  the  top  level  of  the  structure.  The  relation 
between  the  cross-sectional  areas  and  the  moment  of  inertia  are  given  by  I  =  4.62  A2.  Initial 
cross-sectional  areas  for  all  bars  are  30  in2.  Due  to  symmetry  of  the  frame  there  are  only  two 
independent  loading  conditions  in  view  of  the  design  variable  linking. 


Young’s  Modulus 

29.0  x  106  psi 

Poisson’s  Ratio 

0.3 

Weight  Density 

0.283  lb/in3 

Table  5-1:  Material  Properties  for  Steel 
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Attachment  1 

Space  Truss  Finite  Element  Bulk  Data 
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Rectangular  Wing  Finite  Element  Bulk  Data 
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t 

t 

20,  4, 

0.5 

MPC,  100, 

17, 

2, 

-1.0, 

20,  2, 

1.0, 

f 

f 

20,  4, 

-0.5 

MPC,  100, 

18, 

2, 

-1.0, 

20,  2, 

1.0, 

f 

t 

20,  4, 

0.5 

$  Root  connection 

veritically 

MPC,  100, 

13, 

3, 

-1.0, 

20,  3, 

1.0, 

r 

r 

20,  5, 

10.0 

MPC,  100, 

14, 

3, 

-1.0, 

20,  3, 

1.0, 

t 

t 

20,  5, 

10.0 

MPC,  100, 

15, 

3, 

-1.0, 

20,  3, 

1.0 

MPC,  100, 

16, 

3, 

-1.0, 

20,  3, 

1.0 

MPC,  100, 

17, 

3, 

-1.0, 

20,  3, 

1.0, 

r 

t 

20,  5, 

-10.0 

MPC,  100, 

$ 

18, 

3, 

-1.0, 

20,  3, 

1.0, 

f 

t 

20,  5, 

-10.0 

$  Top  Skins 
$ 


CQDMEM1 

10001 

10001 

1 

3 

9 

7 

0.0 

CQDMEM1 

10002 

10001 

3 

5 

11 

9 

0.0 

CQDMEM1 

10003 

10002 

7 

9 

15 

13 

0.0 

CQDMEM1 

10004 

10002 

9 

11 

17 

15 

0.0 
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PQDMEM1 

10001 

2 

0.2000 

PQDMEM1 

e 

10002 

2 

0.2000 

V 

$  Bottom 
e 

Skins 

CQDMEM1 

20001 

20001 

2 

4 

10 

8 

0.0 

CQDMEM1 

20002 

20001 

4 

6 

12 

10 

0.0 

CQDMEM1 

20003 

20002 

8 

10 

16 

14 

0.0 

CQDMEM1 

20004 

20002 

10 

12 

18 

16 

0.0 

PQDMEM1 

20001 

2 

0.2000 

PQDMEM1 

20002 

2 

0.2000 

$ 

$  Spars 


5 

CSHEAR 

30001 

30001 

1 

2 

8 

7 

CSHEAR 

30002 

30001 

3 

4 

10 

9 

CSHEAR 

30003 

30001 

5 

6 

12 

11 

CSHEAR 

30004 

30001 

7 

8 

14 

13 

CSHEAR 

30005 

30001 

9 

10 

16 

15 

CSHEAR 

30006 

30001 

11 

12 

18 

17 

PS HEAR 

c 

30001 

1 

0.0500  0.0 

y 

$  Ribs 

c 

CSHEAR 

40001 

40001 

1 

3 

4 

2 

CSHEAR 

40002 

40001 

3 

5 

6 

4 

CSHEAR 

40003 

40001 

7 

9 

10 

8 

CSHEAR 

40004 

40001 

9 

11 

12 

10 

CSHEAR 

40005 

40001 

13 

15 

16 

14 

CSHEAR 

40006 

40001 

15 

17 

18 

16 

PSHEAR 

e 

40001 

1 

0.0500  0.0 

y 

$  Posts 

e 

y 

CROD 

50001 

50001 

1 

2 

CROD 

50002 

50001 

3 

4 

CROD 

50003 

50001 

5 

6 

CROD 

50004 

50001 

7 

8 

CROD 

50005 

50001 

9 

10 

CROD 

50006 

50001 

11 

12 

CROD 

50007 

50001 

13 

14 

CROD 

50008 

50001 

15 

16 

CROD 

50009 

50001 

17 

18 

PROD 

<! 

50001 

1 

0.01  0. 

00 

0.000 

0.000 

y 

$  Materials 

e 

V 

MAT1 

1 

10.0E+6 

0.30 

0.10 

+MT1 

+MT1  20 

.  0E+3 

15.0E+3 

12  .OE+3 

MAT1 

2 

10.0E+6 

0.30 

0.10 

+MT2 

+MT2  20 

e 

.  OE+3 

15. OE+3 

12. OE+3 

y 

$  Fuselage 

Mass 

V 

CONM2 ,  1,  20,  ,  200.0,  -10.0,  0.0,  0.0 

9  9 

+CONM 

+CONM,  ,  , 

22500. 

,  ,  ,  22500. 

Rectangular  Wing  bulk  data 
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</>  </>  </> 


Weight  to  Mass 
CONVERT,  MASS,  0 


Rectangular  Wing  bulk  data 


Conversion  (for  densities  and  lumped  masses) 
.  00259 
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Cantilevered  Plate  Finite  Element  Bulk  Data 

$ 

$  15  DEGREE  SWEPT,  UNTAPERED  WING.  (SEE  NACA  RM  L55E11  FOR  A 

$  DESCRIPTION  OF  THE  MODEL)  THE  STRUCTURE  IS  MODELLED  USING 

$  15  CQUAD4  ELEMENTS  IN  A  3  X  15  ELEMENT  GRID.  24  GRID  POINTS 

$  ARE  REQUIRED  TO  DEFINE  THE  CORNERS  OF  THE  BENDING  ELEMENTS. 

$ 


$  M= . 45 

EXPERIMENTAL 

RESULTS : 

FLUTTER  VEL  = 

495  EPS  (59 

IN/SEC) 

$ 

c 

FLUTTER  FREQ  = 

120  HZ 

V 

GRID 

1 

0.000 

0 . 000 

0 . 000 

123456 

GRID 

2 

0.690 

0.000 

0 . 000 

123456 

GRID 

3 

1.380 

0.000 

0.000 

123456 

GRID 

4 

2.070 

0 . 000 

0 . 000 

123456 

9 

GRID 

5 

0.296 

1 .  104 

0 . 000 

GRID 

6 

0.986 

1.104 

0.000 

GRID 

7 

1 .  676 

1 . 104 

0.000 

GRID 

c 

8 

2.366 

1 . 104 

0.000 

9 

GRID 

9 

0.592 

2.208 

0.000 

GRID 

10 

1.282 

2.208 

0.000 

GRID 

11 

1 . 972 

2.208 

o.ooc 

GRID 

12 

2.662 

2.208 

0.000 

V 

GRID 

13 

0.888 

3.312 

0.000 

GRID 

14 

1.57  8 

3.312 

0.000 

GRID 

15 

2.268 

3.312 

0.000 

GRID 

<> 

16 

2.958 

3.312 

0.000 

9 

GRID 

17 

1.184 

4.416 

0.000 

GRID 

18 

1 . 874 

4.416 

0.000 

GRID 

19 

2.564 

4.416 

0.000 

GRID 

e 

20 

3.254 

4.416 

0.000 

9 

GRID 

21 

1 .480 

5 . 520 

0.000 

GRID 

22 

2 . 170 

5 . 520 

0.000 

GRID 

23 

2.860 

5.520 

0.000 

GRID 

e 

24 

3.550 

5.520 

0.000 

9 

CQUAD4 

101 

100 

1 

2 

6 

5 

CQUAD4 

102 

100 

2 

3 

7 

6 

CQUAD4 

103 

100 

3 

4 

8 

7 

CQUAD4 

104 

100 

5 

6 

10 

9 

CQUAD4 

105 

100 

6 

7 

11 

10 

CQUAD4 

106 

100 

7 

8 

12 

11 

CQUAD4 

107 

100 

9 

10 

14 

13 

CQUAD4 

108 

100 

10 

11 

15 

14 

CQUAD4 

109 

100 

11 

12 

16 

15 

CQUAD4 

110 

100 

13 

14 

18 

17 

CQUAD  4 

111 

100 

14 

15 

19 

18 

CQUAD4 

112 

100 

15 

16 

20 

19 

CQUAD 4 

113 

100 

17 

18 

22 

21 

CQUAD 4 

114 

100 

18 

19 

23 

22 

CQUAD 4 

115 

100 

19 

20 

24 

23 
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Swept  Wing  Finite  Element  Bulk  Data 


$ 

$  SWEPT  WING  MODEL  FROM 

$  "A  ROOT  LOCUS  BASED  FLUTTER  SYNTHESIS  PROCEDURE"  BY 

$  P.  HAJELA  STANFORD  U. 

$  WITH  A  FLUTTER  CONSTRAINT  AT  SEA  LEVEL  FOR  M=0.80 

$  STRESS  CONSTRAINTS  UNDER  A  4  G  STATIC  AIR  LOAD  AT 

$  25000  FT.  (M  =  1.25)  AND  A  1.5  HZ  LOW.  BOUND  FREQ.  CONSTRNT . 


$ 


GRID 

1 

0.0 

0 . 0 

10.039 

GRID 

2 

0.0 

0.0 

-10.039 

GRID 

3 

72.8345 

0.0 

10.039 

GRID 

4 

72.8345 

0.0 

-10 . 039 

GRID 

5 

145 .6690 

0.0 

10.039 

GRID 

6 

145.6690 

0 . 0 

-10 . 039 

GRID 

7 

53.4758 

116. 667 

9 . 3502 

GRID 

8 

53 . 4758 

116.667 

-9.3502 

GRID 

9 

121 . 1590 

116.667 

9 . 3502 

GRID 

10 

121.1590 

116.667 

-9.3502 

GRID 

11 

188.8430 

116 . 667 

9.3502 

GRID 

12 

188.8430 

116.667 

-9.3502 

GRID 

13 

106.9520 

233.333 

8 .6613 

GRID 

14 

106.9520 

233.333 

-8.6613 

GRID 

15 

169 . 4840 

233.333 

8 . 6613 

GRID 

16 

169.4840 

233.333 

-8 . 6613 

GRID 

17 

232 .0170 

233.333 

8 . 6613 

GRID 

18 

232 .0170 

233.333 

-8 . 6613 

GRID 

1  ; 

160.4280 

350.0 

7 . 9724 

GRID 

20 

160 . 4280 

350.0 

-7 . 9724 

GRID 

21 

217 . 8090 

350.0 

7 . 9724 

GRID 

22 

217 . 8090 

350.0 

-7 . 9724 

GRID 

23 

275.1910 

350.0 

7 . 9724 

GRID 

24 

275.1910 

350 . 0 

-7 . 9724 

GRID 

25 

213 . 9030 

466 . 667 

7.2834 

GRID 

26 

213.9030 

466 . 667 

-7.2834 

GRID 

27 

266.1340 

466 .667 

7.2834 

GRID 

28 

266.1340 

466.667 

-7.2834 

GRID 

29 

318.3650 

466 . 667 

7.2834 

GRID 

30 

318 .3650 

466 .667 

-7.2834 

GRID 

31 

267 . 3780 

583 . 333 

6 . 5945 

GRID 

32 

267.3780 

583.333 

-6 . 5945 

GRID 

33 

314 . 4590 

583 . 333 

6.5945 

GRID 

34 

314 . 4590 

583 . 333 

-6 . 5945 

GRID 

35 

361 . 5390 

583.333 

6 . 5945 

GRID 

36 

361 . 5390 

583 . 333 

-6.5945 

GRID 

37 

320 . 8550 

700 . 0 

5.9055 

GRID 

38 

320 .8550 

700 . 0 

-5 . 9055 

GRID 

39 

362.7840 

700 . 0 

5 . 9055 

GRID 

40 

362.7840 

700 . 0 

-5 . 9055 

GRID 

41 

404 .7130 

700 . 0 

5.9055 

GRID 

42 

404.7130 

700 . 0 

-5.9055 

GRID 

43 

290.7840 

700.0 

0.0 

GRID 

44 

434.7830 

700.0 

0.0 

GRID 

45 

72 .8345 

o 

o 

o 

o 

$ 


( 
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$  BOUNDARY  CONDITION  1 

$ 


MPC ,  101, 

43, 

1, 

O 

1 

37,  1, 

1.0,  ,  MPC4311 

+PC4311, 

38, 

1, 

1.0, 

39,  1, 

1.0,  ,  MPC4312 

+PC4312, 

40, 

1, 

1 . 0 

MPC,  101, 

44, 

1, 

-4.0, 

39,  1, 

1.0,  ,  MPC 4  411 

+PC4411, 

40, 

1, 

1.0, 

41,  1, 

1.0,  ,  MPC4412 

+PC4412, 

42, 

1, 

1 . 0 

MPC,  101, 

43, 

2, 

-4.0, 

37,  2, 

1.0,  ,  MPC4321 

+PC4321, 

38, 

2, 

1.0, 

39,  2, 

1.0,  ,  MPC4322 

+PC4322 , 

40, 

2, 

1.0 

MPC,  101, 

44, 

2, 

-4.0, 

39,  2, 

1.0,  ,  MPC4421 

+PC4421, 

40, 

2, 

1.0, 

41,  2, 

1.0,  ,  MPC4422 

+PC4422 , 

42, 

2, 

1.0 

MPC,  101, 

43, 

3, 

-1 . 

0,  37, 

3,  0.85859,  ,  MPC4331 

+PC4331, 

38, 

3, 

0.85859,  39, 

3,-0.35859,  ,  MPC4332 

+PC4332 , 

40, 

3, 

-0.35859 

MPC,  101, 

44, 

3, 

-1  . 

0,  39, 

3,-0.35859,  ,  MPC4431 

+PC4431, 

40, 

3, 

-0.35859,  41, 

3,  0.85859,  ,  MPC4432 

+PC4432, 

42, 

3, 

0 .85859 

SPC1,  10 

,  123456, 

1, 

THRU, 

6,  45 

SPC1,  10 

f 

456, 

7, 

THRU,  44 

ASET1 ,  100, 

3,  7 

,  9, 

11,  13, 

15,  17, 

ASETA 

+SETA,  19, 

21,  23 

,  25, 

27,  29, 

31,  33, 

ASETB 

+SETB,  35, 

37,  39 

,  41 

$ 

$  BOUNDARY  CONDITION  2 

$ 

MPCADD ,  2101,  101,  201 


MPC, 

201, 

3, 

1, 

1.0, 

45, 

5, 

-10.04 

MPC, 

201, 

3, 

3, 

1.0, 

45, 

3, 

-1.0 

MPC, 

201, 

4, 

1, 

1.0, 

45, 

5, 

10.04 

MPC, 

201, 

4, 

3, 

1.0, 

45, 

3, 

-1.0 

SPC1, 

110, 

1246, 

45 

SPC1, 

110, 

2456, 

1, 

THRU, 

6 

SPC1, 

110, 

456, 

7, 

THRU, 

44 

ASET1 , 

1100, 

3,  7 

,  9, 

11 

,  13, 

15, 

17, 

ASETA 

+SETA, 

19, 

21,  23, 

25, 

27, 

29, 

31, 

33,  ASETB 

+SETB, 

35, 

37,  39, 

41, 

45, 

1,  5 

ASET1,  1100,  5,  45 
SUPORT,  1,  45,  35 

$ 

$  UPPER  AND  LOWER  SKINS  100  -  UPPER,  200  -  LOWER 

$ 


CQDMEM1 

101 

1004 

1 

7 

9 

3 

CQDMEM1 

201 

1004 

2 

8 

10 

4 

CQDMEM1 

102 

1004 

3 

9 

11 

5 

CQDMEM1 

202 

1004 

4 

10 

12 

6 

CQDMEM1 

103 

1004 

7 

13 

15 

9 

CQDMEM1 

203 

1004 

8 

14 

16 

10 

CQDMEM1 

104 

1004 

9 

15 

17 

11 

CQDMEM1 

204 

1004 

10 

16 

18 

12 

CQDMEM1 

105 

1005 

13 

19 

21 

15 

CQDMEM1 

205 

1005 

14 

20 

22 

16 

CQDMEM1 

106 

1005 

15 

21 

23 

17 

CQDMEM1 

206 

1005 

16 

22 

24 

18 

CQDMEM1 

107 

1005 

19 

25 

27 

21 

CQDMEM1 

207 

1005 

20 

26 

28 

22 

Swcpi  Wing  bulk  data 


415 


Attachment  4 


CQDMEMl 

108 

1005 

21 

27 

29 

23 

CQDMEMl 

208 

1005 

22 

28 

30 

24 

CQDMEMl 

109 

1006 

25 

31 

33 

27 

CQDMEMl 

209 

1006 

26 

32 

34 

28 

CQDMEMl 

110 

1006 

27 

33 

35 

29 

CQDMEMl 

210 

1006 

28 

34 

36 

30 

CQDMEMl 

111 

1006 

31 

37 

39 

33 

CQDMEMl 

211 

1006 

32 

38 

40 

34 

CQDMEMl 

112 

1006 

33 

39 

41 

35 

CQDMEMl 

c; 

212 

1006 

34 

40 

42 

36 

MODEL  SUB 

STRUCTURE 

$ 

SHEAR  PANELS:  300 

-  LE, 

35C 

)  -  MID, 

400 

-  TE,  500  -  CHORDWISE 

$ 

AXIAL  RODS 

:  600  - 

INBOARD 

2 

BAYS 

$ 

700  - 

MID  SPAN 

O 

z. 

BAYS 

$ 

c 

800  - 

OUTBOARD 

2 

BAYS 

V 

CSHEAR 

301 

2007 

1 

2 

8 

7 

CSHEAR 

351 

2007 

3 

4 

10 

9 

CSHEAR 

401 

2007 

5 

6 

12 

11 

CSHEAR 

302 

2007 

7 

8 

14 

13 

CSHEAR 

352 

2007 

9 

10 

16 

15 

CSHEAR 

402 

2007 

11 

12 

18 

17 

CSHEAR 

303 

2008 

13 

14 

20 

19 

CSHEAR 

353 

2008 

15 

16 

22 

21 

CSHEAR 

403 

2008 

17 

18 

24 

23 

CSHEAR 

304 

2008 

19 

20 

26 

25 

CSHEAR 

354 

2008 

21 

22 

28 

27 

CSHEAR 

404 

2008 

23 

24 

30 

29 

CSHEAR 

305 

2009 

25 

26 

32 

31 

CSHEAR 

355 

2009 

27 

28 

34 

33 

CSHEAR 

405 

2009 

29 

30 

36 

35 

CSHEAR 

306 

2009 

31 

32 

38 

37 

CSHEAR 

356 

2009 

33 

34 

40 

39 

CSHEAR 

406 

2009 

35 

36 

42 

41 

CSHEAR 

501 

2010 

7 

8 

10 

9 

CSHEAR 

502 

2010 

9 

10 

12 

11 

CSHEAR 

503 

2010 

13 

14 

16 

15 

CSHEAR 

504 

2010 

15 

16 

18 

17 

CSHEAR 

505 

2011 

19 

20 

22 

21 

CSHEAR 

506 

2011 

21 

22 

24 

23 

CSHEAR 

507 

2011 

25 

26 

28 

27 

CSHEAR 

508 

2011 

27 

28 

30 

29 

CSHEAR 

509 

2012 

31 

32 

34 

33 

CSHEAR 

510 

*.012 

33 

34 

36 

35 

CSHEAR 

511 

2012 

37 

38 

40 

39 

CSHEAR 

512 

2012 

39 

40 

42 

41 

CSHEAR 

513 

2010 

1 

2 

4 

3 

CSHEAR 

e 

514 

2010 

3 

4 

6 

5 

•P 

CONROD 

1201 

1 

2 

90 

0.3 

CONROD 

1202 

3 

4 

90 

0.3 

CONROD 

1203 

5 

6 

90 

0.3 

CONROD 

1301 

7 

8 

90 

0.3 

CONROD 

1302 

13 

14 

90 

0.3 

CONROD 

1303 

19 

20 

90 

0 . 3 

CONROD 

1304 

25 

26 

90 

0.3 

CONROD 

1305 

31 

32 

90 

0.3 
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CONROD 

1306 

37 

38 

90 

0.3 

CONROD 

1401 

9 

10 

90 

0.3 

CONROD 

1402 

15 

16 

90 

0.3 

CONROD 

1403 

21 

22 

90 

0.3 

CONROD 

1404 

27 

28 

90 

0.3 

CONROD 

1405 

33 

34 

90 

0.3 

CONROD 

1406 

39 

40 

90 

0.3 

CONROD 

1501 

11 

12 

90 

0.3 

CONROD 

1502 

17 

18 

90 

0.3 

CONROD 

1503 

23 

24 

90 

0.3 

CONROD 

1504 

29 

30 

90 

0.3 

CONROD 

1505 

35 

36 

90 

0.3 

CONROD 

1506 

41 

42 

90 

0.3 

CROD 

601 

6001 

1 

7 

CROD 

602 

6001 

2 

8 

CROD 

603 

6001 

3 

9 

CROD 

604 

6001 

4 

10 

CROD 

605 

6001 

5 

11 

CROD 

606 

6001 

6 

12 

CROD 

607 

6001 

7 

13 

CROD 

608 

6001 

8 

14 

CROD 

609 

6001 

9 

15 

CROD 

610 

6001 

10 

16 

CROD 

611 

6001 

11 

17 

CROD 

612 

6001 

12 

18 

CROD 

701 

7002 

13 

19 

CROD 

702 

7002 

14 

20 

CROD 

703 

7002 

15 

21 

CROD 

704 

7002 

16 

22 

CROD 

705 

7002 

17 

23 

CROD 

706 

7002 

18 

24 

CROD 

707 

7002 

19 

25 

CROD 

708 

7002 

20 

26 

CROD 

709 

7002 

21 

27 

CROD 

710 

7002 

22 

28 

CROD 

711 

7002 

23 

29 

CROD 

712 

7002 

24 

30 

CROD 

801 

8003 

25 

31 

CROD 

802 

8003 

26 

32 

CROD 

803 

8003 

27 

33 

CROD 

804 

8003 

28 

34 

CROD 

805 

8003 

29 

35 

CROD 

806 

8003 

30 

36 

CROD 

807 

8003 

31 

37 

CROD 

808 

8003 

32 

38 

CROD 

809 

8003 

33 

39 

CROD 

810 

8003 

34 

40 

CROD 

811 

8003 

35 

41 

CROD 

e 

812 

8003 

36 

42 

s> 

CONM2 

50001 

7 

20 . 0 

CONM2 

50002 

8 

20.0 

CONM2 

50003 

9 

20.0 

CONM2 

50004 

10 

20.0 

CONM2 

50005 

11 

20.0 

CONM2 

50006 

12 

20.0 

CONM2 

50007 

13 

20.0 

CONM2 

50008 

14 

20.0 

Swept  Wing  bulk  data 
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CONM2 

50009 

15 

20 . 0 

CONM2 

50010 

16 

20.0 

CONM2 

50011 

17 

20 . 0 

CONM2 

50012 

18 

20.0 

CONM2 

50013 

19 

20.0 

CONM2 

50014 

20 

20 . 0 

CONM2 

50015 

21 

20.0 

CONM2 

50016 

22 

20.0 

CONM2 

50017 

23 

20.0 

CONM2 

50018 

24 

20.0 

CONM2 

50019 

25 

20.0 

CONM2 

50020 

26 

20.0 

CONM2 

50021 

27 

20.0 

CONM2 

50022 

28 

20 . 0 

CONM2 

50023 

29 

20.0 

CONM2 

50024 

30 

20.0 

CONM2 

50025 

31 

20.0 

CONM2 

50026 

32 

20 . 0 

CONM2 

50027 

33 

20 . 0 

CONM2 

50028 

34 

20.0 

CONM2 

50029 

35 

20.0 

CONM2 

50030 

36 

20 . 0 

CONM2 

50031 

37 

40 . 0 

CONM2 

50032 

38 

40 . 0 

CONM2 

50033 

39 

40.0 

CONM2 

50034 

40 

40.0 

CONM2 

50035 

41 

40.0 

CONM2 

50036 

42 

40.0 

CONM2 

50037 

43 

40.0 

CONM2 

50038 

44 

40.0 

$ 

$  TRIM  WEIGHT  AT  ROOT  1/4  CHORD  INCLUDING  ROTATIONAL  INERTIA 

$ 

CONM2,  51001,  45,  ,  30000.0,  -36.0,  ,  ,  ,  +CM01 


+CM01, 

$ 

r 

,  3.6E9 

PQDMEM1 , 

1004, 

91, 

0.04 

PQDMEM1 , 

1005, 

91, 

0.04 

PQDMEM1, 

$ 

1006, 

91, 

0.04 

PSHEAR, 

2007, 

90, 

0.04 

PSHEAR, 

2008, 

90, 

0.04 

PSHEAR, 

2009, 

90, 

0.04 

PSHEAR, 

2010, 

90, 

0.04 

PSHEAR, 

2011, 

90, 

0.04 

PSHEAR, 

$ 

2012, 

90, 

0.04 

PROD, 

6001, 

90, 

1 . 0 

PROD, 

7002, 

90, 

1.0 

PROD, 

$ 

8003, 

90, 

1.0 

$  Material 

properties 

$ 

MAT1, 

90, 

10  . 

E6 ,  ,  0.3, 

0.1 

MAT1 , 

91, 

10. 

E6,  ,  0.3, 

0.1,  ,  ,  ,  ABC 

+BC,  60000.0, 

50000 

.0,  30000.0 

$ 

CONVERT, 

MASS, 

2 . 588E 

-3 

Swept  Wing  bulk  dan 
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Attachment  5 

Plane  Frame  Finite  Element  Bulk  Data 


GRDSET, 

t 

,  ,  ,  ,  345 

GRID,  1, 

f 

0  . ,  1200 . 

,  o. 

GRID,  2, 

t 

180  . ,  1200 . 

,  o. 

GRID,  3, 

t 

360  . ,  1200 . 

,  o. 

GRID,  4, 

/ 

0 . ,  1080  . 

,  o. 

GRID,  5, 

r 

180  . ,  1080  . 

,  o. 

GRID,  6, 

f 

360  . ,  1080  . 

,  o. 

GRID,  7, 

r 

0., 960., 

0  . 

GRID,  8, 

/ 

180. ,  960., 

0  . 

GRID,  9, 

/ 

360., 960., 

0  . 

GRID, 10, 

/ 

0  . ,  840  . , 

0  . 

GRID, 11, 

/ 

180  . ,  840  . , 

0  . 

GRID, 12, 

t 

360  . ,  840  . , 

0  . 

GRID, 1 3, 

t 

0  .  ,  720  .  , 

0  . 

GRID, 14, 

t 

180  . ,  720  .  , 

0  . 

GRID, 15, 

! 

360  . ,  720  .  , 

0  . 

GRID, 16, 

r 

0 . ,  600  . , 

0  . 

GRID, 17, 

t 

180 . ,  600  . , 

0  . 

GRID, 18, 

t 

360. ,600., 

0  . 

GRID, 19, 

t 

0  . ,  480  . , 

0  . 

GRID, 20, 

t 

180. ,  480., 

0. 

GRID, 21, 

r 

360 . ,  480  . , 

0. 

GRID, 22, 

t 

0 . ,  360  . , 

0. 

GRID, 23, 

r 

180. ,360., 

0. 

GRID, 24, 

r 

360. ,360., 

0  . 

GRID, 25, 

r 

0  . ,  240  . , 

0  . 

GRID, 26, 

t 

180. ,240., 

0. 

GRID, 27, 

t 

360. ,240., 

0. 

GRID, 28, 

t 

0., 120., 

0. 

GRID, 29, 

r 

180. ,120., 

0. 

GRID, 30, 

r 

360. ,120., 

0  . 

GRID, 31, 

t 

0-,  0., 

0. 

GRID, 32, 

t 

360.,  0., 

0. 

BAROR,  , 

t 

,  ,  0 . ,  0 . 

,  1. 

CBAR,  1, 

i, 

1,  2,  ,  , 

,  1000  . 

CBAR,  2, 

i. 

2,  3,  ,  , 

,  1000. 

CBAR,  3, 

i, 

4  5 

"  /  -J  !  t  t 

,  1000  . 

CBAR,  4, 

i. 

5 ,  6 ,  ,  , 

,  1000  . 

CBAR,  5, 

i. 

7,  8 ,  ,  , 

,1000. 

CBAR,  6, 

i, 

8,  9,  ,  , 

,  1000  . 

CBAR,  7, 

i, 

10,  11,  , 

,  ,1000. 

CBAR,  8, 

i, 

11,  12,  , 

,  ,1000. 

CBAR,  9, 

i. 

13,  14,  , 

,  ,1000. 

CBAR, 10, 

i, 

14,  15,  , 

,  ,1000. 

CBAR, 11, 

i. 

16,  17,  , 

,  ,1000. 

CBAR, 12, 

i. 

t— * 

00 

,  ,1000. 

CBAR, 13, 

i. 

19,  20,  , 

,  ,1000. 

CBAR, 14, 

i. 

20,  21,  , 

,  ,1000. 

CBAR, 15, 

i, 

22,  23,  , 

,  ,1000. 

CBAR, 16, 

i, 

23,  24,  , 

,  ,1000. 

CBAR, 17, 

i, 

25,  26,  , 

,  ,1000. 

CBAR, 18, 

i, 

26,  27,  , 

,  ,1000. 

CBAR, 19, 

i. 

28,  29,  , 

,  ,1000. 

CBAR, 20, 

i, 

29,  30,  , 

,  ,1000. 

CBAR, 21 , 

i, 

1,  4,  , 

, 1000. 

Plane  Frame  bulk  data 
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CBAR, 22, 

1, 

3, 

6,  , 

, 1000  . 

CBAR, 23, 

1, 

4, 

7,  , 

,  1000  . 

CBAR, 24, 

1, 

6, 

9,  , 

,  1000  . 

CBAR, 25, 

1, 

7, 

10,  , 

,  1000  . 

CBAR, 26, 

1, 

9, 

12,  , 

,  1000  . 

CBAR, 27, 

1, 

10, 

13, 

f 

,  ,1000 

CBAR, 28, 

1, 

12, 

15, 

t 

,  ,1000 

CBAR, 29, 

1, 

13, 

16, 

r 

,  ,1000 

’BAR, 30, 

1, 

15, 

CO 

/ 

,  ,1000 

„BAR, 31, 

1, 

16, 

19, 

r 

,  ,1000 

CBAR, 32, 

1, 

CD 

21, 

t 

,  ,1000 

CBAR, 33, 

1, 

19, 

22, 

t 

,  ,1000 

CBAR, 34, 

1, 

21, 

24, 

t 

,  ,1000 

CBAR, 35, 

1, 

22, 

25, 

/ 

,  ,1000 

CBAR, 36, 

1, 

24, 

27, 

/ 

,  ,1000 

CBAR, 37, 

1, 

25, 

28, 

t 

,  ,1000 

CBAR, 38, 

1, 

27, 

30, 

! 

,  ,1000 

CBAR, 39, 

1, 

28, 

31, 

f 

,  ,1000 

CBAR, 40, 

1, 

30, 

32, 

f 

,  ,1000 

FORCE, 

1, 

1, 

/ 

1000  . , 

0 

*  / 

-6 

•  t 

o . 

FORCE, 

1,  2, 

f 

1000  . , 

0 

•  / 

-24 

•  / 

o . 

FORCE, 

1,  3, 

/ 

1000., 

0 

•  / 

-6 

•  r 

0  . 

FORCE, 

1,  4, 

t 

1000., 

0 

•  / 

-12 

*  ! 

0. 

FORCE, 

1,  5, 

t 

1000  . , 

0 

•  r 

-48 

*  r 

0. 

FORCE, 

1, 

6, 

r 

1000., 

0 

•  1 

-12 

•  f 

0  . 

FORCE, 

1,  7, 

r 

1000., 

0 

•  r 

-12 

•  t 

0  . 

FORCE, 

1,  8, 

t 

1000  . , 

0 

*  t 

-48 

*  r 

0. 

FORCE, 

1, 

9, 

r 

1000., 

0 

•  / 

-12 

•  r 

0. 

FORCE, 

1,  io, 

f 

1000  . , 

0 

•  t 

-12 

•  r 

0. 

FORCE, 

1,  11, 

i 

1000  . , 

0 

•  r 

-48 

•  / 

0  . 

FORCE, 

1,12, 

t 

1000., 

0 

•  r 

-12 

•  / 

0. 

FORCE, 

1,13, 

f 

1000., 

0 

•  r 

-12 

•  r 

0. 

FORCE, 

1,  14, 

t 

1000  . , 

0 

•  r 

-48 

•  r 

0  . 

FORCE, 

1,  15, 

t 

1000., 

0 

•  r 

-12 

*  f 

0. 

FORCE, 

1,16, 

r 

1000., 

0 

•  t 

-12 

•  f 

0. 

FORCE, 

1,  17, 

t 

1000  . , 

0 

•  f 

-48 

'  r 

0. 

FORCE, 

1,  18, 

r 

1000., 

0 

•  f 

-12 

•  r 

0  . 

FORCE, 

1,  19, 

r 

1000., 

0 

•  f 

-12 

•  t 

0  . 

FORCE, 

1,20, 

f 

1000., 

0 

•  r 

-48 

•  t 

0  . 

FORCE, 

1,21, 

f 

1000  . , 

0 

•  r 

-12 

•  / 

0. 

FORCE, 

1,22, 

r 

1000., 

0 

•  r 

-12 

■  t 

0  . 

FORCE, 

1,23, 

r 

1000  . , 

0 

•  i 

-48 

*  t 

0  . 

FORCE, 

1,24, 

/ 

1000  . , 

0 

•  r 

-12 

•  t 

0  . 

FORCE, 

1,25, 

f 

1000., 

0 

•  r 

-12 

•  t 

0. 

FORCE, 

1,26, 

t 

1000., 

0 

•  f 

-48 

•  r 

0  . 

FORCE, 

1,27, 

r 

1000., 

0 

•  r 

-12 

•  t 

0  . 

FORCE, 

1,28, 

f 

1000  . , 

0 

•  / 

-12 

•  r 

0  . 

FORCE, 

1,29, 

f 

1000., 

0 

•  t 

-48 

•  r 

0  . 

FORCE, 

1,30, 

/ 

1000., 

0 

•  r 

-12 

•  r 

0  . 

FORCE, 

2,  1, 

t 

1000., 

0 

•  r 

-4  . 

5, 

0  . 

FORCE, 

2,  2, 

t 

1000., 

0 

•  r 

-18 

•  r 

0  . 

FORCE, 

2, 

3, 

t 

1000  . , 

0 

•  t 

-4  . 

5, 

0  . 

FORCE, 

2, 

4, 

r 

1000  . , 

0 

•  t 

-9 

*  f 

0  . 

FORCE, 

2,  5, 

t 

1000., 

0 

•  t 

-36 

•  r 

0  . 

FORCE, 

2,  6, 

r 

1000., 

0 

•  t 

-9 

•  t 

0  . 

FORCE, 

2,  7, 

t 

1000., 

0 

•  f 

-9 

•  t 

0  . 

FORCE, 

2,  8, 

r 

1000  . , 

0 

•  r 

-36 

•  t 

0  . 

FORCE, 

2,  9, 

r 

1000., 

0 

•  t 

-9 

•  t 

0  . 

Mane  F  rame  bulk  dam 
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f 


FORCE, 

2,  10, 

9 

1000  . 

/ 

0., 

-9  . 

0  . 

FORCE, 

2,  11, 

9 

1000  . 

9 

o.. 

-36. 

0  . 

FORCE, 

2,  12, 

9 

1000  . 

9 

0., 

-9  . 

0  . 

FORCE, 

2,  13, 

9 

1000  . 

9 

0., 

-9  . 

0  . 

FORCE, 

2,  14, 

1 

1000  . 

9 

0., 

-36  . 

0. 

FORCE, 

2,  15, 

9 

1000  . 

9 

0., 

-9. 

0  . 

FORCE, 

2,16, 

9 

1000. 

9 

0., 

-9. 

0. 

FORCE, 

2,  17, 

9 

1000  . 

9 

0., 

-36  . 

0  . 

FORCE, 

2,18, 

t 

1000. 

9 

0., 

-9. 

0  . 

FORCE, 

2,19, 

t 

1000  . 

9 

0., 

-9. 

0  . 

FORCE, 

2,20, 

9 

1000  . 

9 

0., 

-36. 

0. 

FORCE, 

2,21, 

t 

1000  . 

9 

0., 

-9. 

0. 

FORCE, 

2,22, 

9 

1000. 

9 

0., 

-9  . 

0  . 

FORCE, 

2,23, 

9 

1000. 

9 

0., 

-36. 

0  . 

FORCE, 

2,24, 

t 

1000  . 

t 

0., 

-9  . 

0  . 

FORCE, 

2,25, 

9 

1000. 

9 

0., 

-9. 

0  . 

FORCE, 

2,26, 

r 

1000  . 

9 

0., 

-36. 

0  . 

FORCE, 

2,27, 

9 

1000  . 

9 

0., 

-9  . 

0  . 

FORCE, 

2,28, 

/ 

1  O00  . 

9 

0., 

-9. 

0. 

FORCE, 

2,29, 

9 

1000  . 

9 

0., 

-36. 

0  . 

FORCE, 

2,30, 

9 

1000. 

9 

0., 

-9. 

0  . 

FORCE, 

2,  1, 

9 

1000  . 

9 

3.18 

,  o. 

0  . 

FORCE, 

2,  4, 

9 

1000  . 

9 

6.23 

,  o. 

0. 

FORCE, 

2,  7, 

9 

1000. 

9 

6.08 

,  0. 

0  . 

FORCE, 

2,10, 

9 

1000  . 

9 

5.92 

,  o. 

0. 

FORCE, 

2,13, 

9 

1000  . 

9 

5.74 

,  o. 

0  . 

FORCE, 

2,16, 

9 

1000. 

9 

5.54 

,  o. 

0. 

FORCE, 

2,19, 

9 

1000  . 

9 

5.30 

,  o. 

0. 

FORCE, 

2,22, 

9 

1000  . 

9 

5.00 

,  o. 

0. 

FORCE, 

2,25, 

9 

1000. 

9 

4.61 

,  0. 

0. 

FORCE, 

2,28, 

9 

1000. 

9 

4.00 

,  0. 

0. 

FORCE, 

3,  1, 

9 

1000. 

9 

o., 

-4.5 

0. 

FORCE, 

3,  2, 

9 

1000. 

9 

o.. 

-18. 

0  . 

FORCE, 

3,  3, 

9 

1000  . 

9 

0., 

-4.5 

0  . 

FORCE, 

3,  4, 

9 

1000  . 

9 

o.. 

-9  . 

0  . 

FORCE, 

3,  5, 

9 

1000. 

9 

0., 

-36. 

0  . 

FORCE, 

3,  6, 

9 

1000  . 

9 

o.. 

-9. 

0  . 

FORCE, 

3,  7, 

9 

1000  . 

9 

0., 

_9  _ 

0  . 

FORCE, 

3,  8, 

9 

1000. 

9 

0., 

-36. 

0. 

FORCE , 

3,  9, 

9 

1000  . 

9 

o.. 

-9. 

0. 

FORCE, 

3,  10, 

9 

1000. 

9 

o.. 

-9  . 

0  . 

FORCE, 

3,  11, 

9 

1000  . 

9 

0., 

-36. 

0  . 

FORCE, 

3,  12, 

9 

1000  . 

9 

0., 

-9. 

0. 

FORCE, 

3,  13, 

9 

1000. 

9 

0., 

-9. 

0. 

FORCE, 

3,  14, 

9 

1000  . 

9 

0., 

-36. 

0  . 

FORCE, 

3,  15, 

9 

1000  . 

9 

0-, 

-9. 

0  . 

FORCE, 

3,16, 

9 

1000  . 

9 

0., 

-9. 

0  . 

FORCE, 

3,17, 

9 

1000  . 

9 

0., 

-36. 

0. 

FORCE, 

3,18, 

9 

1000  . 

9 

0., 

-9. 

0  . 

FORCE, 

3,  19, 

9 

1000  . 

9 

0., 

-9. 

0. 

FORCE, 

3,20, 

9 

1000  . 

9 

0., 

-36. 

0. 

FORCE, 

3,21, 

9 

1000. 

9 

0., 

-9. 

0  . 

FORCE, 

3,22, 

9 

1000  . 

9 

0., 

-9. 

0. 

FORCE, 

3,  23, 

9 

1000  . 

9 

0., 

-36. 

0. 

FORCE, 

3,24, 

9 

1  000  . 

9 

0., 

-9  . 

0  . 

FORCE, 

3,25, 

9 

1000  . 

9 

0., 

-9. 

0  . 

FORCE, 

3,26, 

9 

1000  . 

9 

0., 

-36  . 

0. 

FORCE, 

3,27, 

9 

1000  . 

9 

0., 

-9  . 

0  . 

Plane  Frame  bulk  data 
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FORCE,  3,28, 
FORCE,  3,29, 
FORCE,  3,30, 
FORCE,  3,  3, 

FORCE,  3,  6, 

FORCE,  3,  9, 

FORCE,  3,12, 
FORCE,  3,  15, 
FORCE,  3,18, 
FORCE,  3,21, 
FORCE,  3,24, 
FORCE,  3,27, 
FORCE,  3,30, 
SPC1,  6,  126 

MAT1 ,  1,  29. 

(MAT  I ,  2  4.  (  3 


Plane  Frame  bulk  dala 


,  1000.,  o.,  -9., 

0  . 

,  1000.,  0.,  -36., 

0  . 

,  1000.,  0.,  -9., 

0  . 

,  1000.,  -3.18,  0., 

0. 

,  1000.,  -6.23,  0., 

0. 

,  1000.,  -6.08,  0., 

0  . 

,  1000.,  -5.92,  0., 

0  . 

,  1000.,  -5.74,  0., 

0  . 

,  1000.,  -5.54,  0., 

0  . 

,  1000.,  -5.30,  0., 

0. 

,  1000.,  -5.00,  0., 

0. 

,  1000.,  -4.61,  0., 

0  . 

,  1000.,  -4.00,  0., 

0. 

,  31,  32 

+6,  ,  0.3,  0.283,  , 

,  ,  +MAT1 

,  24.(3,  24.(0 
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